Resonant Inelastic X-ray Scattering on Spin-Orbit Coupled Insulating Iridates 
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We determine how the elementary excitations of iridium-oxide materials, which are dominated 
by a strong relativistic spin-orbit coupling, appear in Resonant Inelastic X-ray Scattering (RIXS). 
Whereas the RIXS spectral weight at the L2 x-ray edge vanishes, we find it to be strong at the 
Z,3-edge. Applying this to S^IrCU, we observe that RIXS, besides being sensitive to local doublet 
to quartet transitions, meticulously maps out the strongly dispersive delocalized excitations of the 
low-lying spin-orbit doublets. 
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Materials containing ions with an orbital degeneracy 
show a plethora of physical effects related to the sponta- 
neous lifting of this degeneracy pQ . The two degeneracy- 
breaking mechanisms in Mott insulators traditionally 
considered are the cooperative Jahn- Teller effect [2], in- 
volving a change in lattice symmetry, and the superex- 
change interactions [3 , intertwining long-range ordering 
of orbital and magnetic degrees of freedom. A third 
and less explored possibility exists in systems containing 
heavy ions, where strong relativistic spin-orbit coupling 
dominates the spin-orbital physics [4]. 

The strong spin-orbit interaction can cause entirely 
new kinds of ordering that are of topological nature. This 
was recently proposed for certain iridium-oxides [6J, 
members of a large family of iridium-based materials. 
Na2lr03, for instance, is predicted to be a topologi- 
cal insulator exhibiting the quantum spin Hall effect at 
room temperature [5j • The topologically non-trivial state 
arises from the presence of complex hopping integrals, 
resulting from the unquenched iridium orbital moment. 
This system can also be described in terms of a Mott 
insulator, with interactions between the effective irid- 
ium spin-orbital degrees of freedom goverened by the 
Kitaev-Heisenberg model (7HH] • In the pyrochlore iridates 
A2IX2O7 (where A is a 3+ ion), a quantum phase tran- 
sition from a topological band insulator to a topological 
Mott insulator has been proposed as a function of the 
electron-electron interaction strength [5J [TOl E] . 

To establish whether and how such novel phases are 
realized in iridium oxides it is essential to probe and 
understand their spin-orbital ordering and related ele- 
mentary excitations. In this context it is advantageous 
to consider the structurally less complicated, single-layer 
iridium perovskite Sr 2 Ir04. This material is in many 
respects the analog of the high-T c cuprate parent com- 
pound La2Cu04 [8 . Structurally it is identical, with the 
obvious difference that the Ir 5d valence electrons are, as 
opposed to the Cu 3d electrons, very strongly spin-orbit 
coupled. The similarity cuts deeper, however, as the low 
energy sector of the iridates is spanned by local spin-orbit 



doublets with an effective spin of 1/2, which reside on a 
square lattice and interact via superexchange - a close 
analogy with the undoped cuprates. This observation 
motivates doping studies of Sr2lr04 searching for super- 
conductivity [12 US] • Experimentally, however, far less 
is known about the microscopic ordering and excitations 
in iridates than in cuprates. Inelastic neutron scatter- 
ing, which can in principle reveal such properties, is not 
possible because Ir is a strong neutron absorber [T3] and, 
moreover, crystals presently available tend to be tiny. As 
a consequence not even the interaction strength between 
the effective spins in the simplest iridium-oxides is estab- 
lished: estimates for Sr 2 Ir04, for instance, range from 
~50 meV [H] to ~110 meV Q2]. 

In this Letter we show that while for iridates neutron 
scattering falls short, photon-in photon-out scattering in 
the form of resonant inelastic x-ray scattering (RIXS) [T5] 
fills the void: RIXS at the iridium L-edge offers direct ac- 
cess to the excitation spectrum across the Brillouin zone, 
enabling one to measure the dispersion of elementary 
magnetic excitations. Besides the low energy magnons 
related to long-range order of the doublets, RIXS will 
also reveal the dynamics of higher energy, doublct-to- 
quartet, spin-orbit excitations. This allows to directly 
test theoretical models for the excitation spectra and ex- 
tract accurate values of the superexchange and spin-orbit 
coupling constants J and A, respectively. 

7r 4+ ionic ground state. — In the iridium-oxides Ir + 
ions are located in octahedra of oxygen ions, splitting the 
5d levels by ~ 3 eV into e g and t 2g orbitals [T7] . Because 
this crystal field splitting is an order of magnitude larger 
than the spin-orbit coupling A, the t 2g levels do not hy- 
bridize much with the e g orbitals, and a t\ q configuration 
is established [T5] . The symmetry of the t 2g ground state 
is in principle governed by three factors: superexchange 
interactions, additional lattice-induced crystal field split- 
tings, and relativistic spin-orbit coupling 19J. The su- 
perexchange J [8] is estimated to be about an order of 
magnitude smaller than A w 0.4 eV '20 . An elongation 
of the octahedra along the z axis [3T] favors a ground 
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state where the hole is in the xy orbital. But experimen- 
tal data strongly favor the spin-orbit coupling scenario 
over the lattice splitting scenario, however (4J [18] . 

The orbital degree of freedom of the hole can be de- 
scribed by an effective angular momentum 1 = 1, related 
to the true orbital angular momentum by 1 = — L |22j . 
The orbital eigenstates of l z are described by the anni- 
hilation operators rfo,±ij defined in terms of the real t 2g 
wave functions by the relations d yz — — {d\ — d_i)/V2, 
d zx = i{d\ + d-i)/y/2, d xy = do- When the spin-orbit 
coupling term is projected to the t 2g subspace, it be- 
comes — Al • S. Tetragonal lattice distortions can also be 
included, and the Hamiltonian for a single Ir ion is [8] 
H = — Al • S — A/ 2 , with A > for elongation along 
the z axis. The six eigenstates group into three Kramers 
doublets, described by the fermions /, g and h with an- 
nihilation operators 
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and with energies w/ = A/(v2 tan#), ui g = —A — A/2 and 
uj h = -(Atan0)/V2, where tan 29 = 2\/2A/(A - 2A). 
For A = 0, which corresponds to the cubic, isotropic 
situation sin# = y/l/3 and cos 9 = y/2/3. For A/A <C 1, 
the hole's ground state doublet is {\xy t) , \xy 4-)}- For 
A/ A 3> 1, the eigenstates are characterized by the total 
effective angular momentum J e g = 1 + S. In the ground 
state, the hole occupies the / doublet (J e g = 1/2), which 
is separated by an energy of 3 A/2 from the J e g = 3/2 
quartet, which splits into the g and h doublets. 

The remaining, two-fold ground state degeneracy can- 
not be removed by Jahn- Teller distortions because the 
two states in the ground state Kramers doublet have ex- 
actly the same charge distribution. Superexchange cou- 
pling, however, is present in all iridates, and couples the 
local doublets, thus dictating the low energy collective 
dynamics of the material. 

RIXS cross section. — Resonant Inelastic X-ray Scat- 
tering is particularly suited to probe higher energy mag- 
netic excitations and dispersions, as demonstrated in the 
cuprates 23-26 . The recent increase in brilliance of 
the new generation synchrotron X-ray sources allows for 
spectra with tremendous resolving powers, enough to re- 
solve dispersion up to a few tens of meV. In RIXS [IB|l2"7]. 
a photon with momentum Kk, energy and polariza- 
tion e is scattered to Hk', av and e', losing momentum 
hq = hk—hk 1 and energy cu = lu^—lo^ to the sample, cjk 
is tuned to a certain atomic resonance of the material un- 
der study, greatly enhancing the scattering cross section. 
In our case, that will be the Ir L-edge: the 2p core elec- 
tron is excited into the empty 5d t 2g state. After a very 
short time, another electron from the t 2g levels can fall 
back to the core hole under the emission of an outgoing 



X-ray. The system is left in an excited state, whose en- 
ergy and momentum are taken from the scattered X-ray 
photon, which is measured. 

The RIXS cross section is described by the Kramers- 
Heisenberg equation [35], where the photon absorption 
and subsequent emission are governed by the dipole op- 
erator V = J2r e lk ' R r • e acting on all electrons of an 
Ir 4+ ion at site R. The phonon polarization is e. 

The intermediate state has a filled shell (5d t 2g ), so the 
dominant multiplet effect comes from the core orbital's 
spin-orbit coupling A: the 2p core states split into J = 
1/2 (the L 2 edge) and J = 3/2 states (the L 3 edge). Since 
the L 2 and the L 3 edge are separated by 1.6 keV @], their 
interference is negligible, given the much smaller lifetime 
broadening of a few eV [29]. Because the 2p core states 
have the same angular momenta as the 5d t 2g states, we 
can describe them with the three fermions F, G and H , 
analogous to Eq. Q, where we replace (d yz , d zx , d xy ) by 
{Vx,Py,Pz) and the parameters A, A and 9 by A, S and 0. 
The tetragonal distortion S is expected to be very small 
for the deep 2p core states. 

The lifetime broadening at the Ir L edge is still 
quite large compared to the dynamics of the 5d elec- 
trons [18l [29] - Therefore, we make the fast collision ap- 
proximation Ei + ftw-k — E n + iT iT 30 . The sum 
over n in the Kramers-Heisenberg equation can be per- 
formed, and comprises the core states of either the L 2 
or the L3 edge. In second quantization, the dipole op- 
erators are r • e = Y, a ,p, a ( 5c U r \ 2 Pp) ' e d l* Pp a + h - c - 
which can be denoted as (D2 + -D3) +h.c, where D 23 are 
the local dipole transition operators for the L 2 and L 3 
edge, respectively. The RIXS amplitude becomes A q oc 
(/I E R e iq ' R [D t (e / *)D(e)] R |0), where R runs over all Ir 
sites and the RIXS intensity 7 q (w) = J2f |Aj| 2 <5(w — Ef)- 
We rewrite the inelastic scattering operator as 

&D = £ [Bliftf.+BiiflU + BlUlg. 
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Integrating out the core hole degree of freedom one ob- 
tains at the L 2 edge the intra-doublet scattering matrix 
elements = -sin 2 (6» - 6)4* e ff and B f J s = 0. The 
doublet-quartet excitation matrix elements of the spin- 
orbit multiplet are = sin(6» - 6) cos 6 e'*e s , B f a % = 
-(-l) CT sin(6> - 9) sin9 e'*e„, = -\{-l) a sin2(6> - 
<d)e'st s and b£ = 0, where (-1) CT is 1 for a = t and 
— 1 for it = J,. Further, = e + and = e_, with 
e± = (e x ± ie y )/^/2. In the case of dominant spin-orbit 
coupling, 9 = 0. Since all matrix elements at the L 2 
edge are proportional to sin(# — 6), the inelastic scatter- 
ing intensity vanishes completely in this case, in addition 
to a vanishing of the elastic intensity [4]. 

At the L3 edge, however, RIXS is fully allowed. The 
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FIG. 1: (Color online.) RIXS spectra of Sr 2 Ir0 4 at the Ir L 3 
t2 g edge. The left panel shows the spectrum for incoming n 
polarization, and the right one for incoming a polarization. 
The outgoing polarization is not measured. The intra-doublet 
excitations are broadened by J/10. 
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FIG. 2: Vertical cuts through Fig. [T] The left panel shows 
spectra at several q for incoming n polarization, and the right 
one for incoming a polarization. 
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6)e£e ff - cos 2 (9 e'*e z and B S J S = \{-tf sin26» (e£e 
e'*e a ) for the intra-doublet ones, and Bl 9 — cos(( 
8) sin 6 e' z *e s , b{% = (-1) CT cos(0-9) cos 9 e'*e s , B& 
\{-lY[sm2{6 - 9)4% - sin20 (e' CT *e CT - e' z *e z ) s , - 
— sin 2 9 e' z * e a — cos 2 8 e'£ e z for the doublet-quartet exci- 
tations. For excitations within the J c g — 1/2 doublet, 
the scattering operator can be rewritten in terms of the 
effective angular momentum operator, which in the limit 
A/ A -c 1 takes the particularly simple form D\d z — 
| (e'* • e 1 + P • J cff ), where P x = i{e'*e z - e'* e y ) and its 
cyclic permutations P VjZ are polarization factors. Here, 
the first term corresponds to elastic scattering while the 
P • Jeff term gives rise to inelastic scattering. 

RIXS on Sr 2 IrOi. - Up to this point, the dis- 
cussion is general and applies to all materials with an 
Ir 4+ ion in an octahedral crystal field. Calculation of 
the RIXS spectra for a particular iridate is straightfor- 
ward given the Hamilonian that captures the interac- 
tions between the Ir degrees of freedom. In Sr2lr04 
the effective low energy Hamiltonian is obtained from 
the spin-orbital superexchange for the triply degener- 
ate t 2g orbitals (Eq. (3.11) from Ref. [T!5]) by project- 
ing it on the low energy Kramers doublet. In the case of 
strong spin-orbit coupling, one finds a Heisenberg Hamil- 
tonian for these pseudo-spin- 1/2 states, with weak dipo- 
lar anisotropy due to Hund's rule coupling. The rotation 
of the octahedra around the z axis over an angle a ks 11° 
introduces a Dzyaloshinsky-Moriya interaction, but after 
an appropriate spin rotation the Hamiltonian remains of 
Heisenberg type [8]. 

At the Ir L 3 edge, excitations within the J c ff = 1/2 
doublet can be described in terms of Holstein-Primakoff 
bosons. The single- and double-magnon intensities are, 
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with Uk = (l/yr^ 2 + i) 1 / 2 /^, Vk = (i/yr^_ 

l) 1 / 2 sign(7 k )/v / 2 and 7 k = (cosfe^ +cosfc y )/2. A re- 
markable difference with L-edge RIXS on cuprates [23] 
is that the large spin canting, reflected in the appreciable 
value of a, causes the presence of spectral weight in the 
center of the Brillouin-zone, at q = 0. 

Transitions from J c ff = 1/2 to 3/2, which are at an 
energy of |A, are expected to show a less pronounced 
q dependence. The crystal field splitting of the quar- 
tet states is probably too small to resolve with current 
RIXS instruments, so we give the integrated intensity of 
all these excitations: I {9+h) oc 2 + \e' ■ e| 2 - \e'* ■ e\ 2 . 
The polarization terms cancel unless both incoming and 
outgoing X-rays are circularly polarized. 

Computed RIXS spectra. — We now evaluate the dif- 
ferent contributions to the RIXS intensity. Single and 
double magnon contributions i^ 1 ' 2 ) and those from the 
JcS — 1/2 to 3/2 excitations I^ 9+h ^> are presented for the 
specific case of a 90° scattering angle with the scatter- 
ing plane perpendicular to the Ir0 2 layers, and q in the 
first (2D) Brillouin zone. The resulting cross sections are 
shown in Fig s.|l]andg where we used |A/J = 10. The 
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FIG. 3: (Color online.) Spectral weight of different excita- 
tions. The units on the vertical axis are chosen such that the 
Jeff = 1/2 to 3/2 excitation has a spectral weight of unity. 



low-energy intra doublet excitations with A J ff = show 
a very distinct magnon dispersion, the intensity of which 
is strongly varying with q. The doublet-quartet transi- 
tions with locally AJ e g = 1 are at |A, corresponding 
to 0.5 - 0.6 eV [20 . This implies they match in energy 
the large spectral weight charge modes observed in opti- 
cal absorption in Sr2lr04 [3T]. Even if the local multi- 
plet excitations are not optically active themselves, there 
will be strong mixing of the spin-orbit excited state with 
inter-site charge excitations across the Mott gap. This 
causes a derealization of the doublet-quartet mode on 
the scale of the intersite hopping t. The dispersion and 
momentum dependent spectral weight modulations that 
this causes is beyond the present model calculations; here 
it only reflects in the use of an effective broadening of the 
doublet-quartet mode with t, corresponding to about 4 J. 

To summarize, we have determined the effective scat- 
tering operators for direct RIXS at the L-edge of Ir 4+ 
ions in an octahedral crystal field. In the physical limit 
of strong spin-orbit coupling, the RIXS spectral weight 
at the L 2 vanishes, but it is strong at the L 3 -edge. Ap- 
plying this to Sr2lrC>4, we find that RIXS can map out 
the strongly dispersive single- and double-magnon excita- 
tions of the low-lying doublet and is in addition sensitive 
to the doublet-quartet excitations at an energy of §A, 
which strongly mix with delocalized charge modes. This 
shows that RIXS can accurately determine the material 
parameters A and J of iridates and is an excellent tool 
to probe their low-energy elementary excitations, testing 
and characterizing proposed long-range order and topo- 
logical phases. 

We thank H. Takagi, J. H. Kim, B. J. Kim and M. van 
Veenendaal for valuable discussions. 
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